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Abstract 

The current paper is aimed in getting more insight on three main 
points concerning large-scale astrophysical systems, namely: (i) for- 
mulation of tensor virial equations from the standpoint of analytical 
mechanics; (ii) investigation on the role of systematic and random 
motions with respect to virial equilibrium configurations; (iii) extent 
to which systematic and random motions are equivalent in flatten- 
ing or elongating the shape of a mass distribution. The tensor virial 
equations are formulated regardless from the nature of the system and 
its constituents, by generalizing and extending a procedure used for 
the scalar virial equations, in presence of discrete subunits (Landau 
& Lifchitz 1966). In particular, the self potential-energy tensor is 
shown to be symmetric with respect to the exchange of the indices, 
{Epot)pg = {Epot)qp- Then the results are extended to continuous 
mass distributions. The role of systematic and random motions in 
collisionless, ideal, self-gravitating fluids, is analysed in detail includ- 
ing radial and tangential velocity dispersion on the equatorial plane, 
and the related mean angular velocity, fi, is conceived as a figure ro- 
tation. R3 fluids are defined as ideal, self-gravitating fluids in virial 
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equilibrium, with systematic rotation around a principal axis of iner- 
tia. The related virial equations are written in terms of the moment of 
inertia tensor, Ipq, the self potential-energy tensor, (-Epot)pg, and the 
generalized anisotropy tensor, C,pq (Caimmi & Marmo 2005; Caimmi 
2006a) . Additional effort is devoted to the investigation of the prop- 
erties of axisymmetric and triaxial configurations. A unified theory of 
systematic and random motions is developed for R3 fluids, taking into 
consideration imaginary rotation (Caimmi 1996b, 2006a). The effect 
of random motion excess is shown to be equivalent to an additional 
real or imaginary rotation, respectively, inducing flattening (along the 
equatorial plane) or elongation (along the rotation axis). Then it is 
realized that a R3 fluid always admits an adjoint configuration with 
isotropic random velocity distribution. In addition, further constraints 
are established on the amount of random velocity anisotropy along the 
principal axes, for triaxial configurations. A necessary condition is for- 
mulated for the occurrence of bifurcation points from axisymmetric to 
triaxial configurations in virial equilibrium, which is independent of 
the anisotropy parameters. A particularization is made to the spe- 
cial case of homeoidally striated Jacobi ellipsoid, and some previously 
known results (Caimmi 2006a) are reproduced. 

keywords - scalar virial theorem; tensor virial theorem; collisionless, 
self- gravitating ideal fluids; random velocity distributions; stellar sys- 
tems. 

1 Introduction 

Large-scale celestial objects, such as stellar systems, galaxy clusters, and 
(non baryonic) dark matter haloes predicted by current ACDM cosmologies, 
may safely be represented as collisionless, ideal self-gravitating fluids. The 
related flow equation takes the same formal expression as in their coUisional 
counterpart, with the exception that the pressure force is generalized in terms 
of a stress tensor, allowing different rms velocities along different directions 
[e.g., Binney & Tremaine 1987 (hereafter quoted as BT87), Chap. 4, §2]. 
Accordingly, collisionless fluids can be flattened equally well by rotation (with 
respect to a selected axis) and/or anisotropic random velocity distribution 
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i.e. anisotropic pressure (e.g., Caimmi 2006a, hereafter quoted as C06^). 
In fact, giant elliptical galaxies exhibit a negligible amount of (systematic) 
rotation, and their shape is mainly due to anisotropic pressure (e.g., Bertola 
& Capaccioli 1975; lUingworth 1977, 1981; Schechter & Gunn 1979; BT87, 
Chap. 4, §36). 

Collisionless fluids of astrophysical and cosmological interest range over 
about ten decades in mass, from globular clusters to galaxy clusters say, pro- 
vided gaseous i.e. coUisional component may safely be neglected. Therefore, 
it seems necessary to investigate the role of systematic and random motions 
in making virializcd collisionless fluids. To this respect, the virial theorem 
in tensor form may be a useful tool. According to the standard procedure, 
the tensor virial equations are determined along the following steps (e.g., 
Binney 1978, 2005; Wiegandt 1982a,b; BT87, Chap. 4, §3): (a) start with 
the collisionless Boltzmann equations; (b) derive a set of moment equations; 
(c) integrate the above set of moment equations, under some simplifying 
assunptions. 

To the (limited) knowledge of the author, no attempt can be found in the 
literature where (i) the tensor virial equations are formulated from the stand- 
point of analytical mechanics; (ii) the role of systematic and random motions 
is clearly stated, and (iii) the equivalence between systematic and random 
motions in flattening or elongating the boundary, is clearly established. 

Concerning (i), the tensor virial equations could be determined regard- 
less from the nature of the system and its constituents, by generalizing and 
extending a procedure used for the scalar virial equations [Landau & Lifchitz 
1966 (hereafter quoted as LL66), Chap. II, § 10]. With regard to (ii), prelimi- 
nary considerations reported in previous attempts [Caimmi 1996a,b; Caimmi 
& Marmo 2005 (hereafter quoted as CM05); C06] should be further improved 
and developed. In dealing with (iii), the definition of imaginary rotation al- 
lows an interpretation of rms velocity excess in terms of systematic rotation 
around a fixed principal axis (Caimmi 1996b; C06; Caimmi 2006b^) which 
could be inserted in the context under discussion. A detailed investigation 
on the above mentioned points makes the aim of the present attempt. 

The current paper is organized as follows. A general formulation of the 

more extended file including an earlier version of the above quoted paper is available 
at the arxiv electronic site, as astro-ph/0507314. 

more extended file including an earlier version of the above quoted paper is available 
at the arxiv electronic site, as astro-ph/0507314. 



3 



tensor virial theorem which holds, in particular, for the gravitational inter- 
action, is provided in Sect. 2. The properties of R3 fluids, defined as ideal, 
self-gravitating fluids in virial equilibrium, rotating around a principal axis, 
are studied in Sect. 3, and a unified theory of systematic and random motions 
is provided in Sect. 4, where some general relations are particularized to the 
special case of homeoidally striated Jacobi ellipsoids, and previously known 
results (Caimmi 1996a,b; C06) are reproduced. Some concluding remarks 
are drawn in Sect. 5, and a few arguments are treated with more detail in 
the Appendix. 



2 The tensor virial theorem 

A general procedure used for the formulation of the scalar virial theorem 
(LL66, Chap. II, § 10) shall be followed here in the derivation of the tensor 
virial theorem. Let us take into consideration a mechanical system made of 
N particles, referred to an inertial frame. Let {xi)r, {vi)r, be the position 
and velocity components related to i particle, and rrii the mass, 1 < i < N , 
1 < r < 3. 

The kinetic-energy tensor: 

1 ^ 

{Ei^n)pq ^ -J2mi{Vi)p{Vi)g ; (1) 

^ i=i 

is a function of 2N or variables, (t'i),., 1 < < A, 1 < p < 3, 1 < g < 3, 
for selected p and q, according if the tensor components are diagonal or non 
diagonal, respectively. The kinetic-energy tensor is manifestly symmetric 
with respect to the indices: 

(-£'kin)p(j — (-E'kin)qp j (2) 

and the trace is the kinetic energy: 

3 N 3 

s=l ^ i=l s=l 

which is a function of 3A^ variables, {vi)r, l<i<N,l<r<3. The first 
partial derivatives are: 

(Pi)r = ^ ^i('^i)r (4) 
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where {pi)r, 1 < i < iV, 1 < r < 3, is the impulse component of i particle 
(e.g., LL66, Chap. II, §7). 

The combination of Eqs. (1) and (4) yields: 



N 



(5) 



i=l 



which is equivalent to: 



2(-E'kin)pg — 



N 



ijq 



U=l 



N 



(6) 



or, using Newton's equations (e.g., LL66, Cap. I, §5): 



2(-E'kin)pq' — 



N 



Li=l 



N 



+ E( 

i=l 



X. 



pot 



ijp 



(7) 



ijq 



where Epot[{xi)r], l<i<N,l<r<3,is the self potential energy. 

The last term on the right-hand side of Eq. (7) defines a tensor, the trace 
of which is usually named the virial of the system (Clausius 1870). In the 
author's opinion, it would be better to quote the virial and its parent tensor as 
the virial potential energy and the virial potential-energy tensor, respectively. 

If the self potential energy is a homogeneous function of the coordinates, 
of degree x, then the following relation holds: 



which, in turn, implies: 

(-£'pot)pg[C('^j)r] — C''''(-^pot)pq'[(^j)r] j 

where {Epot)pq is defined as the self potential-energy tensor: 



TV 



{Epot) pq 



dE, 



pot 



(8) 



(9) 



(10) 



and C is a generic real number, provided C{xi)ri 1<'?<A^, l<'r<3, 
belonges to the domain of -Spot- 
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With regard to the self potential energy, the Euler theorem reads: 



3 N dEpot 
and the combination of Eqs. (10) and (11) yields: 

3 



(11) 



(-S'pot )ss — Epot ; 



(12) 



as expected. 

The substitution of Eq. (10) into (7) yields: 



_d 
dt 



N 



.1=1 



— 2(£^kin)pqr — x{Epot)pq '■> 



(13) 



and the sum of Eq. (13) with its counterpart where the indices, p and q, are 
interchanged, reads: 



N 



-nJ2ii^i)piPi)l + i^i)qiPi)p\ = 2[(^kin)pg + (^kin)gp]-X[(^pot)p5 + (^pot)gp] • 
^'^ i=l 

(14) 

Let us define the moment of inertia tensor^ [e-g-, Chandrasekhar 1969 
(hereafter quoted as C69), Chap. 2, §9; BT87, Chap. 4, §3]: 



N 



— ^ j '^iiS^i)'piS^i)q 1 
i=l 



(15a) 
(15b) 



s=l 



where / is the total moment of inertia of the system, with respect to the 
centre of mass. Owing to Eq. (4), the first temporal derivative is: 



'■pi 



= [(^Op(Pi)<z + i.^i)qipi)p] ; 



(16) 



^In this formulation, the moment of inertia with respect to a coordinate 3iXlS ^ ^ IS 
Ir = Ipp + Iqq, r ^ p ^ q. For a different formulation where Ir = Irr, r = 1,2, 3, see LL66 
(Chap. VI, §32). 
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and the combination of Eqs. (2), (14), and (16), yields: 

Ipg = 4(£'kin)pQ - x[(-E'pot)pg + {Epot)qp\ ', (17) 

on the other hand, the difference of Eq. (13) with its counterpart where the 
indices, p and q, are interchanged, reads: 

d ^ 

-TT [i^i)piPi)<l - i^i)qiPi)p] = -X[{Epot)pq - {Epot)qp] ■ (18) 

i=l 

owing to Eq. (2). 

With regard to the vectors, fi[{xi)i, {xi)2, {xi)3] and Pi[{pi)i, {pi)2, (p^a], 

— * 

and to the vector product, Ji — fi x p^, the sum on the left-hand side of 
Eq. (18) reads (e.g., Spiegel 1968, Chap. 2.2, §§11-12): 

N N 

J2 [(^i)p(Pi)q - i^i)qiPi)p] = ' (^ X Pi) 

1=1 1=1 
N ^ 

— ver^{xr) -^^Ji — ver^(xr) • J — Jr ; (19) 

i=l 

where ver^(x,.) is the versor, or unit vector, parallel to the coordinate axis, 
Xr, r ^ p ^ q, and J is the total angular moment of the system. 
The combination of Eqs. (18) and (19) yields: 

^ = -x[(^pot)p. - (^pot) J ; (20) 

and the conservation of angular momentum, which always holds for isolated 
systems (e.g., LL66, Chap. 2, § 9), implies the symmetry of the self potential- 
energy tensor with respect to the exchange of the indices: 

{Epot)pq = (-E'pot)gp ; (21) 
and Eq. (17) takes the form: 

-^ipq — '^{Ekin)pq — X{Epot)pq (22) 

which makes the virial equations of the second order (for the special case of 
gravitational interaction, x — —1, see e.g., C69, Chap. 2, § 11; BT87, Chap. 4, 
§3). 
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The further restriction: 

ipq^O ] 1<P<3; l<g<3; (23) 

makes Eqs. (22) reduce to: 

2iE^n)pq - xiEpot)pq ^ ; l<p<3; 1<?<3; (24) 

which is the expression of the virial theorem in tensor form^. Strictly speak- 
ing, it holds when the moment of inertia tensor has a linear dependence 
on time, Ipg — kpgt, where kpg are constants. The special case, kpg — 0, 
l<j9<3,l<g<3, is related to dynamical or hydrostatic equilibrium 
(e.g., BT87, Chap. 4, §3). 

An alternative restriction is that the first time derivatives of the moment 
of inertia tensor are bounded, as: 

\ipg{t)\ <Mpg ■ l<p<3; l<g<3; (25) 

where Mpg are convenient real numbers. Accordingly, it can be seen that the 
time average of the second time derivetives of the moment of inertia tensor 
are null (e.g., LL66, Chap. II, §10): 

4g = 0; l<p<3; l<g<3; (26) 
which makes Eqs. (22) reduce to: 

2(^kin)p,-x(^pot)p, = ; l<p<3; l<g<3; (27) 

where time averages are calculated over a sufficiently long (ideally infinite) 
period (e.g., LL66, Chap. 2, §10). In presence of periodic motions (e.g., a 
homogeneous sphere undergoing coherent oscillations), time averages can be 
calculated over a single (or a multiple) period. 

For sake of simplicity, in the following the tensor virial theorem shall 
be expressed by Eqs. (24) where the kinetic-energy and self potential-energy 

'*Somc authors prefer a more general formulation, expressed by Eqs. (22) (e.g., BT87, 
Chap. 4, §3). On the other hand, a more restricted formulation, expressed by Eqs. (24), 
has a closer connection with the scalar virial theorem, which explains the choice adopted 
here. 
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tensors are to be intended as instantaneous or time averaged, according if 
the restriction defined by Eq. (23) or (26) holds. 

The particularization of Eqs. (24) to diagonal components, after summa- 
tion on both sides, produces: 



which is the expression of the virial theorem in scalar form (e.g., LL66, 
Chap. II, § 10). Special cases are (a) Newtonian and Coulombian interaction, 
X = — 1, and (b) Hookeian interaction, x = 2. If the system is in dynamical 
or hydrostatic equilibrium, mean values coincide with instantaneous values. 

The above results are quite general and hold regardless from the nature 
of the system and its constituents, provided no dissipation and/or external 
interaction occur. With regard to a specified system, the sole restrictions 
to be made are (i) the evolution takes place within a finite region of the 
phase hyperspace i.e. < |(ic.j),.| < M^., < |(^-^i)r| < M^^^ 1 < < A^, 
1 < r < 3, where M^., M^., are convenient real numbers, and (ii) the self 
potential energy is a homogeneous function of the 3iV coordinates, of degree 



In dealing with continuous matter distributions instead of mass points, 
the particle mass, m.^, has to be replaced by the mass within an infinitesimal 
volume element, dm = p(a;i, a;2, a^s) dxi da;2 dxs, where p is the density, and 
an integration has to be performed over the whole volume, instead of a sum- 
mation on the coordinates related to all the particles. For further details see 
e.g.. Limber (1959). Accordingly, the kinetic-energy tensor and the kinetic 
energy attain their usual expressions (C69, Chap. 2, §9): 



on the other hand, the self potential-energy tensor and the self potential 
energy read: 



2-Bkin — X-Epot — ; 



(28) 




(29) 



(30) 




(31) 



(32) 
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where is a charge density, V is a potential function, defined as the tidal 
potential energy acting on the unit charge (related to the interaction), placed 
at the point under consideration. For further details, see Appendix A. 
Let us define the total-energy tensor, as: 

Epq = (-E'kin)pg + (-^pot)pg j (^3) 

owing to Eqs. (3) and (12), the related trace: 

3 

E = ^ Ess — E]^in + Epoi ; (34) 

s=l 

is the total energy. 

The combination of Eqs. (24), (33), and (28), (34), respectively, yields: 

iEkin)pq = ^ Epq ] (35) 

2 

{Epot)pq = I o^PI ; (36) 

A, ' 



for tensor components, and: 



-Ekin — 77^-^ ; i^'^) 

2 

Epot = rrr^^ ' ^^^^ 



for tensor traces. 



3 Systematic and random motions 

3.1 Basic ideas 

Let a collisionless, self-gravitating fiuid be referred to an inertial frame, 
(Oa;ia;2a;3), where (without loss of generality) the origin coincides with the 
centre of mass. The number of particles within an infinitesimal hypervolume 
of the phase hyperspace at the time, t, is: 

d'^Af = f{xi,X2,X3,Vi,V2,V3,t)dxidx2dx3dvidv2dv3 ; (39) 
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where / > is the distribution function. The number of particles within an 
infinitesimal volume of the ordinary space at the time, t, is: 

dxidx2dx3 j j j f{xi,X2,X3,vi,V2,V3,t)dvidv2dv3 ; (40) 

where the integration has to be performed over the whole volume in veloc- 
ity space. The number density related to the infinitesimal volume element, 
d^S — dxi dx2 dxs, at the time, t, is: 

d^M f f f 

n{Xi,X2,X3,t) = = J J J f{Xi,X2,X3,Vi,V2,V3,t)dVidV2dV3 ; (41) 

if, in addition, the total particle number, J\f, and the total mass, M, are 
conserved, then the following normalization conditions hold: 

///// /■/'(^i'^2,a;3,t;i,W2,V3,t)da;ida;2dx3dt;idi;2di;3 = A/" ; (42) 
p{xi ,X2,X3,t) dxi dx2 dx3 = M ; (43) 



where p is the mass density of the infinitesimal volume element, d^S, and 
the integrations have to be carried over the whole hypervolume in phase 
hypcrspace and the whole volume in ordinary space, respectively. 

From a physical point of view, the volume element is finite instead of 
infinitesimal, but still containing a large amount of particles which, on the 
other hand, is negligible with respect to the total number. Accordingly, the 
following relations hold: 

I <^ AAf{xi,X2,X3,t) <^N' ; max{mi) <^ AM{xi,X2,X3,t) <^ M ; 

(44) 

where AM and AM represent the particle total number and total mass within 
the volume element, d^S", at the time, t, and rrii is the mass of i particle, 
1 < i < AJ\f. The related total mass, AM, may be expressed as: 

AM 

AM{xi,X2,X3,t) = Y,mi ; (45) 
1=1 

and the mean particle mass within the volume element, d^^S", at the time, t, 
reads: 

_/ ,^ AM(a;i,a;2,a;3,t) , . 

'^^''^''^''^'^= AAfix„X2,X3,t) ' 
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according to the general definition of aritlimetic mean. 

From tlie standpoint of a continuous mass distribution, tlie following 
changes have to be made: AA/'(a;i, X2, X3, — > d^jV'; AM{xi,X2,Xs,t) — > 
p{xi,X2, X3, t) d^S; and Eq. (46) takes the form: 

_/ .X / p{xi,X2,X3,t) 

m{x„ X2, xs, t) = p{x„ X2, xs, t)^^ ^^^^^^^^ ; (47) 

in terms of mass density and number density. 

With regard to Eqs. (42) and (43), equivalent expressions are: 

/ / /^(^i'^2,a^3,*)da:^id3:^2da:^3 =^ ; (48) 
/// '^*^^i'^2,a:^3,^) da;ida;2 dx^ 

= J J J m{xi,X2,X3,t)n{xi,X2,X3,t)dxidx2dx3 ^ M ; (49) 

and the division of both sides of Eq. (49) by their counterparts in Eq. (48), 
yields: 

M , , 

m=j^ ; (50) 

where, owing to the theorem of the mean, m is the particle mass averaged over 
the whole volume. Total mass and particle number conservation imply a time 
independent mean particle mass, m. If, in addition, particles with different 
masses are uniformly distributed throughout the whole volume, then the 
mean particle mass within a generic volume element, d^S", equals the mean 
particle mass within the boundary, as: 

m{xx,X2,X3,t) =m ; (51) 

and the system may be considered, in any respect, as made of M identical 
particles of mass m. Accordingly, Eq. (47) reads: 

p(xi,X2,X3,i) =m n(xi,X2,X3,i) ; (52) 

which implies direct proportionality between mass density and number den- 
sity. 
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Using again the theorem of the mean, let us define the mean velocity 
component, v^, and the mean product velocity component, within a 

generic infinitesimal volume element, d^^S" = dxidx2dx3, at the time, t, as: 

_/ .X J J J f{xi,X2,X3,Vi,V2,V3,t)VpdVidV2dV3 

Vp{xi,X2,X3,t) ^ ——— — ; 53 

j j J j{xi,X2,X3,Vi,V2,V3,t) dt;i dt;2 d^s 

/ ,x J J J f{xi,X2,X3,Vi,V2,V3,t)VpVgdVidV2dV3 

VpVq{xi,X2,X3,t) ^ — , — - — ; (54) 

J J J f{Xi,X2,X3,Vi,V2,V3,t) dvi dv2 dvs 

or, using Eq. (41): 

,^ J J J f(xi,X2,X3,Vi,V2,V3,t)VpdVidV2dV3 
Vp{x„X2,X3,t)- n{x,,X2,X3,t) ' ^^^^ 

. _ 1 1 1 f{Xl,X2,X3,Vi,V2,V3,t)VpVqdVidV2dV3 

VpVg[Xi,X2,X3,t) - ^ ^ A ' 

n[xi, X2, X3, t) 

in terms of the number density, n. 

Let us define the distribution function in the velocity space: 

,x f(xi,X2,X3,Vi,V2,V3,t) 

F{xi,X2,X3,Vi,V2,V3,t) = ; (57) 

n{Xi,X2,X3,t) 

which, owing to Eq. (41), satisfies the normalization condition: 

J J j F{xi,X2,X3,vi,V2,V3,t)dvidv2dv3 ^ 1 ; (58) 
and the substitution of Eq. (57) into (55) and (56) yields: 

V^{xi,X2,X3,t) = J J J F{xi,X2,X3,Vi,V2,V3,t)VpdVidV2dV3 ; (59) 
V^{xi,X2,X3,t) ^ j j j F{xi,X2,X3,Vi,V2,V3,t)VpVqdVidV2dV3 ; (60) 

in terms of the distribution function, F. 

Prom a statistical standpoint, the distribution function, F, may be inter- 
preted as a probability density in the velocity space, where F{xi, X2, X3, vi, V2, 
V3, t) dvi dv2 dv3 represents the probability of finding a particle inside the vol- 
ume element, d^S, at the time, t, with velocity components in the range, 
Vr ^ dVr/2, r = 1,2,3. In this view, the velocity components, Vp, and the 
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product velocity components, VpVq, may be considered as random variables. 
According to the general definition of variance and covariance (e.g., Oliva & 
Terrasi 1976, Chap. II, §2.7), the following relations hold: 

{vlr = {v;? + crl ; (61) 
{vpVgY = v;v* + ; (62) 

where the asterisk denotes the expectation value of the related distribution, 
al is the mathematical (intended as opposite to empirical) variance and 
is the mathematical covariance. According if two random variables, a and b, 
are independent, correlated, or anticorrelated, the mathematical covariance, 
aab, is null, positive, or negative, respectively. 

With regard to an infinitesimal volume element, d'^S*, at the time, t, the 
first term on the right-hand side of Eqs. (61) and (62) is related to the velocity 
components of the centre of mass, while the second term is related to velocity 
components with respect to the centre of mass. 

Expectation values and mathematical variances and covariances are a 
priori quantities which cannot be determined by data collections. The related 
observables are arithmetic means and empirical variances and covariances 
(e.g., Oliva & Terrasi 1976, Chap. IV, §4.3), and Eqs. (61) and (62) translate 
into: 

R) = (^)' + < ; (63) 

Vp^ = % + '^vpvg ; (64) 

where the notation of variances and covariances has been left unchanged, 
for sake of simplicity. Finally, Eq. (63) may also be conglobed into Eq. (64), 
using the definitions: 

where cr^^y^^ may be considered as the generic component of an empirical 
covariance tensor. If velocity components arc independent, cr^upVq = Spq^vp^ 
where Spg is the Kronecker symbol, and Eqs. (63) and (64) merge into: 

= + Spga^p ; (66) 

where a^^. = (Ty^^^ = a^^ for simplifying the notation and considering arr 
as velocity dispersions related to random motions, with regard to a generic 
infinitesimal volume element, d^^*, at the time, t. 
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Having the centre of mass of the system been chosen as origin of the 
(inertial) reference frame, mean (over the whole volume) velocity components 
along coordinate axes are necessarily null: u;: — 0; (e.g., LL66, Chap. II, § 8). 
Accordingly, a^^ = (v^), 0"^^ — 0, p q, and the random square velocity 
tensor, (7^^ = VpV^, is diagonal. 

3.2 Radial and tangential velocity dispersion on the 
equatorial plane 

Let us define the axes of the system as the intersection between the related 
volume and principal axes of inertia, and semiaxes the distances between 
the top axes and the centre of mass. In general, semiaxes on opposite sides 
are different in length. Let the coordinate axis, x^, containing the axis, 03, 
be chosen as rotation axis, without loss of generality, as the reference frame 
may arbitrarily be oriented, provided the origin coincides with the centre of 
mass. Let us define the principal plane, (0x1X2), as the equatorial plane of 
the system. 

With regard to a particle placed at a position, r = {xi,X2, X3), and moving 
at a velocity, v — {vi,V2,V3), let feq — {xi,X2) and Veq — {vi,V2,) be the 
related projections onto the equatorial plane, see Fig. 1. 

Cartesian velocity components may be expressed as the algebraic sum of 
radial and tangential velocity projection on the related direction, as: 

vi ^ VeqCos{a + (j)) = {vw)i + {v^)i ^ Vu,cos(f) - v^s'iikJ) ] (67a) 
V2^VeqSm{a + (f)) ^ {vy,)2 + {v^)2^VwSm4) + v^cos(f) ; (67b) 

conversely, radial and tangential velocity components may be expressed as 
the algebraic sum of cartesian velocity projections on the related direction, 
as: 

Vu, ^ Veqcosa ^ {vi)u, + {v2)w ^ viCOS(f) + V2sm(l) ; (68a) 
■u<^ = ■Ueq sin a = (■ui)0 + (■U2)</, = -Vising + ■^2 cos ; (68b) 

where {Vfj)r = [v • ver^(/i)] ver^(//) • ver^{xr); {vr)^ = [v ■ yef^{xr)]yef^{xr) ■ 
ver^(//); /i — w,(f>; r — 1,2; and ver^(c?) is the unit vector with positive 
orientation, along the d direction. 

With regard to a generic infinitesimal volume element, d^S, at the time, 
t, radial and tangential velocity components, defined by Eqs. (68), may be 
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X2 




Figure 1: 



OP = re, = {xl + a;i)i/2 = w ; 
RQ = v^ = sin a ; 
SQ = (v^)i = —V4, sin 4> ; 
PT = = —VujCos4> ; 

= fPR ; a = RPQ ; 
UQ = Vi = t)eq cos(a + <t>) ; 



PR = = COS a ; 
RS = {v^)2 = V4, cos rj> ; 
TR = {v^)^ = v„ sin ^ ; 
a + <t> = P"QU ; 
PU = V2 = Keq sin(a + 0) ; 



Figure 1: Radial {vyj) and tangential (f^) velocity components on the equa- 
torial plane, (Oa;ia;2). The projection onto the equatorial plane of the vector 
radius, r, and the velocity, is denkfted as Teq and Veq, respectively, with 
regard to a generic particle. 



considered as random variables. Owing to a theorem of statistics^, the ex- 
pectation values of the related distributions read: 



— vl cos + ^2 sin ; 
v^ — —vl sin (f) + V2 cos (j) ; 



(69a) 
(69b) 



similarly, the expectation values of the distributions depending on the ran- 



and using the general definitions expressed by Eqs. (61) and (62), the related 
mathematical variances read: 



= i^lT ~ i^w)^ = ^li + si^^ + '^^viv2 COS sin ; (71a) 

= {'"IT - {^lY = ^l^ sin^ + al^ cos^ - 2(7^,^^,2 sin cos ; (71b) 



where the mathematical covariance, cTv-^wz^ is null provided the related velocity 
components, v\ and ^2, are independent. The validity of the relations: 



can easily be checked. 

In terms of the related observables, arithmetic means and empirical vari- 
ances and covariances, Eqs. (69), (70), and (71) translate into: 



^Let mi, 7722, m„, be random variables and /i(mi)dmi, /2(m2)dTO2, 
fni^n) drUn, related distributions, m = J2k=i ^knik an additional random variable, where 
ak are coefficients, and f{m)dm a related distribution. Then the expectation value, m*, 
is expressible via the above linear combination of the expectation values, mj, m* , 

sic \ * 

as: m = Efe=iafc"i^- 





(72) 

(73) 
(74) 



Vw = Vl COS + f 2 sin : 
v^ = — ^^sin0 + IJ^cos0 




(75b) 
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i^w) = ("^i) cos^ + {V2) sin^ 6 + 2viV2 cos sin ; (76a) 
(vj,) — (vf) sin^ + (^1) cos^ — 2lJ]l^sin cos ; (76b) 



= (^^) - ('"wf = 0-^1 cos^ + (7^2 sin^ + 2(7^1^2 cos sin ; (77a) 

= i'^l) - {"f^y = sin^ + (7^2 cos^ - 2(T^j„2 sin cos ; (77b) 

where the notation of variances and covariances has been left unchanged, for 
sake of simphcity. The vahdity of the relations: 

{v^r + {v^r = {vlf + {v^Y ; (78) 



K) + {vl) = {v!) + {vl) ; (79) 



< + < = < + < ; (80) 

can easily be checked. 

If velocity components are independent, (Ty^vq = ^pqO'v ' -^is. (77) reduce 

to: 

'^ww — ^11 COS + (T22 sin ; (81a) 

where a^r = cr^^y^ = a^^ for simplifying the notation and considering (Trr 
as velocity dispersions related to random motions with regard to a generic 
infinitesimal volume element, d^S", at the time, t. 

With regard to the whole volume, S*, at the time, t, let us define positive 
and negative equatorial radial velocity components, directed outwards 

and inwards, respectively, and positive and negative equatorial tangential 
velocity components, v^f,, as related to counterclockwise and clockwise motion, 
respectively, around the rotation axis. 

Owing to the above mentioned theorem of statistics, the following rela- 
tions hold for for expectation values and mathematical variances related to 
the distributions depending on radial and tangential velocity components on 
the equatorial plane: 

'"w^J^J J J[Vw{Xl,X2,X3,t)]*p{xi,X2,X3,t)dXidX2dX3 ] (82a) 

^^^Jll I /b</'(^i'^2,a;3,t)]*p(a;i,a;2,a;3,t)dxidx2dx3 ; (82b) 
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('"wT = J^J J J[vl{Xi,X2,X3,t)]*p{xi,X2,X3,t)dXidX2dX3 ] (83a) 
,2\* _ ^ f f f L.2 



^<py = I I I [v4,{xi,X2,X3,t)]*p{xi,X2,X3,t)dxidx2dx3 ; (83b) 



< = (vlr - Kf ; (84a) 

< = (^J)* - (^;)' ; (84b) 

where the vahdity of Eqs. (72), (73), and (74) is left unchanged. 

In terms of the related observables, arithmetic means and empirical vari- 
ances, Eqs. (82), (83), and (84), translate into: 

^^=-^yy j v;^{xi,X2,X3,t)p{xi,X2.,X3,t)dxidx2dx3 ; (85a) 

"^^J^JJ Jv^{Xl,X2,X3,t)p{Xi,X2,X3,t)dXidX2dX3 ] (85b) 

'"^^Jll I J ^'U^^^^'^^^3,t)p{xi,X2,X3,t)dxidx2dx3 ; (86a) 

vl{xi,X2,X3,t)]p{xi,X2,X3,t)dXidX2dX3 ] (86b) 



vl 



1 

M 



(^L = = i^l) - {vwf ; (87a) 

2 _ 2 _ 7~T^ 



{V4>) ; (87b) 



where the validity of Eqs. (78), (79), and (80) is left unchanged. 

With regard to the angular velocity, Q, and the related moment of inertia, 
I3, the counterparts of Eqs. (85b) and (86b) read: 

= — y y y f^(xi,X2,X3,i)w^p(xi,a;2,X3,i)dxida:2da;3 ; (88) 

Y J J J ^^i^i^^2^^3,t)w'^p{xi,X2,X3,t)dxidx2dx3 ; (89) 

aj,^W-n' ; (90) 
j j j w'^p{xi,X2,X3,t)dxidx2dx3 ; (91) 
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where the angular velocity, ^l, may be conceived as a figure rotation, in 
the sense that the mean is null when performed in a reference frame in 
rigid rotation at the same rate. To this respect, particles with different 
masses must uniformly be distributed within the region of phase hyperspace 
accessible to the system. 

Owing to Eqs. (84b), (89), and (90), the following relation holds: 



where the mean angular velocity, fl, and the empirical variance, un, are 
related to systematic and random motions, respectively, around the rotation 
axis. Accordingly, Eq. (92) may be sphtted as: 



which express the contribution of systematic and random motions along the 
equatorial plane, in terms of tangential and angular velocities. 

The mean radial velocity component, u^, is related to the motion of the 
centre of mass along the equatorial plane. On the other hand, the centre of 
mass coincides with the origin of the coordinates, which implies — and. 



The mean tangential velocity component, v^, and the empirical variance, 
are related to systematic and random motions, respectively, around the 
rotation axis. The diagonal components of the kinetic-energy tensor may be 
expressed in terms of the above mentioned contributions, as: 



(92) 



(93a) 
(93b) 




T'kk ^\I^sys) kk ~l~ (-^rdm)A;fe i 



(94) 



where k = w^cj) in the case under discussion, and: 




(95a) 
(95b) 




(96b) 



(96a) 
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where the indices, sys and rdm, denote systematic and random motions, 
respectively. 

The global contribution: 

r^^ = ^W ; (97) 

depends only on the mass distribution, via the moment of inertia, I^, related 
to the rotation axis, x^, and the tangential velocity component on the equa- 
torial plane, via the rms angular velocity, Q^, regardless from the amount of 
systematic and random motions along the direction under discussion (e.g., 
Meza 2002; C06). 

The contribution of the kinetic-energy tensor component, T^^, to the 
kinetic-energy tensor components, Tn and T22, owing to Eqs. (96) and (97), 
is: 

{T<t>4>)<ig^ll<i<P^ ; g = l,2 ; (98) 
1 —2 

[{T,ys),p<t>\qq = -Iqq^ ] q = l,2 ; (99) 

m^^UU^l^i'ii^-'^^) ; ? = 1,2 ; (100) 

/.-///w(-i,-.,-3,t)d.,d..d.3 ; (101) 
/3 = /ii + /22 ; (102) 

where Ipg is the moment of inertia tensor. 

With regard to equatorial radial kinetic-energy tensor components, the 
combination of Eqs. (80), (90), and (93) yields: 

^L = ^n + ^i2-^(n^-n') ; (103) 

and the contribution of the kinetic-cncrgy tensor component, T^w, to the 
kinetic-energy tensor components, Tu, T22, owing to Eqs. (94) and (102), is: 



(T^«,),, = ^Ma^^ - l/,,(02 - O') ; q=l,2 ; (104) 



to be used together with Eq. (98). 
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It is worth noticing that anisotropic random velocity components distri- 
butions, aTT 7^ ctrr, are not necessarily related to the shape of the system, 
while aww 20-33 are (e.g., BT87, Chap. 4, § 3), where T, R, and W, denote 
tangential, radial, and equatorial velocity components, respectively. In fact, 
a spherically symmetric mass distribution could, in principle, allow purely 
radial or circular orbits. 

3.3 Virial equilibrium configurations 

The particularization of Eq. (22) to Newtonian interaction, x = — 1, after 
combination with Eqs. (98) and (104), allows the formulation of the virial 
equations for the case under discussion. The result is: 

^4 = 7,,n' + M(7j, + (£;pot),, ; ? = 1,2 ; (105a) 

^/33 = Ma|3 + (£;pot)33 ; (105b) 

where IJ3 = 0, the system centre of mass having been chosen as origin of the 
reference frame. 

The virial equations of the second order, expressed by Eqs. (22), in par- 
ticular Eqs. (105), imply the validity of the following assumptions. 

(i) The mechanical system under consideration (in particular a coUisionless, 

self-gravitating fluid) is isolated (e.g., LL66, Chap. I, § 5), which implies 
angular moment conservation (e.g., LL66, Chap. II, §9). 

(ii) The potential energy is a homogeneous function of the coordinates with 
degree, x (iii particular, % = — 1). 

The validity of the further assumption: 

(iii-a) The generic component of the moment of inertia tensor depends lin- 
early on time, according to Eq. (23); 

or, alternatively: 

(iii-b) The first time derivative of the generic component of the moment of 
inertia tensor is a bounded function, according to Eq. (25); 
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makes Eqs. (105) reduce to: 



7,,Q +Malg + {Ep,t),g^O ; g = l,2 ; (106a) 
Mal^ + {Epot)33^0 ; (106b) 

where the variables are to be intended as instantaneous or averaged over a 
sufficiently long time, according if assumption (iii-a) or (iii-b), respectively, 
has been chosen. 

A more general formulation of Eqs. (106), which includes instantaneuos 
configurations under assumption (iii-b), is: 

7,,n' + MC,X + (^pot),, = ; g = l,2 ; (107a) 
MC33<7' + (^pot)33 = ; (107b) 
^af^ + ai^ + Tls ; (107c) 

C.. = ^%^ = % ; r = 1,2,3 ; (107d) 



Cii + C22 + C33 = t^ = ^ = C ; (i07e) 



;;-2 



where Crr be conceived as anisotropy parameters (CM05, C06), E'rdm is 
the random kinetic energy, and Erdm is the effective random kinetic energy 
i.e. the right amount needed for an instantaneous configuration to satisfy 
Eqs. (106). Generalized anisotropy parameters lower or larger than ^/3 imply, 
respectively, lack or excess of random motions along the related direction. On 
the other hand, the ratios: 

C.. = t^ = ^ ; r = 1,2,3 ; (108a) 

Cii + C22 + C33 = 1 ; (108b) 

may be conceived as effective anisotropy parameters (CMOS, C06). 

The parameter, ^, may be conceived as a virial index, where C = 1 cor- 
responds to null virial excess, 2A£'].dm = 2(i^rdm — -^'rdm), which does not 
necessarily imply a relaxed configuration^, C > 1 to positive virial excess, 

^For instance, a homogeneous sphere undergoing coherent oscillations exhibits C > 1 

at expansion turnover and ^ < 1 at contraction turnover. Then it necessarily exists a 
configuration where C = 1 which, on the other hand, is unrelaxed. 
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and C < 1 to negative virial excess. The special case, C = 1) niakes Eqs. (107) 
reduce to (106). 

For sake of simplicity, let us define ideal, self-gravitating fiuids, rotating 
around an axis, 03, for which Eqs. (107) hold i.e. assumptions (i), (ii), and 
(iii) above are valid, as R3 fluids. 

The combination of Eqs. (107a) and (107b) yields: 

7,,=n'-^(Epot)33 + (^pot),g = ; g = l,2 ; (109) 

S,33 

to get further insight, let us express the self potential-energy tensor, (£'pot)rr) 
and the moment of inertia tensor, I^r, in terms of dimensionless tensors, Vrr 
and Trr, respectively, as: 

{Epot)rr = Vrr (^pot) = V 1,2,3 ;(110) 

a Q, 
Igg = Ma%g h = M tt^X^ ^=1,2 ; (111) 



.2W ' 

in addition, let us define the rotation parameter: 

where p — M/ S is the mean density of the system. In the special case of solid- 
body rotation (f2 = f2), Eq. (113) reduces to a notation used for polytropes 
(e.g., Jeans 1929, Chap. IX, §232; Chandrasekhar & Leboviz 1962), and in 
the limit of ellipsoidal homogeneous configurations (p = p), Eqs. (113) re- 
duces to a notation used for MacLaurin spheroids and Jacobi ellipsoids (e.g., 
Jeans 1929, Chap. VIII, §§189-193; C69, Chap. 5, §32, Chap. 6, §39). 

Owing to Eqs. (110), (111) and (112), Eq. (109) may be formulated in 
terms of dimensionless parameters, as: 

(Caan. - C.T'aa) - <332',, = ; ? = 1,2 ; (114) 

which admits real solutions provided the inequality: 

; g = l,2 ; (115) 

saa / aa 



24 



is satisfied, that is the natural extension to R3 fluids of its counterparts 
related to axisymmetric, relaxed mass distributions (Wiegandtl982a,b) and 
homeoidally striated ellipsoids (C06). Imaginary solutions correspond to 
imaginary rotation parameters i.e. imaginary rotation. 
The combination of Eqs. (114) yields after some algebra: 

2^22(C33^11 - (11^33) = ^ll(C33^22 " (22^^33) ; (116) 

or equivalently: 

C33(2'22Pll - 2:11^22) = ^33(22261 - 2:11(22) ; (117) 

which make alternative expressions of the constraint related to virial equilib- 
rium. 

3.4 Axisymmetric and triaxial configurations 

An explicit expression of the rotation parameter, v, can be derived from 
Eqs. (114), as: 

^^C33^^^_C^ ; g = l,2 ; (118) 
which, in turn, allows an explicit expression of anisotropy parameter ratios, 

Cpp/Cqqi 

? = 1,2 ; (119) 

(120) 

and the combination of Eqs. (107e) and (119) yields: 

C33 _ '^33 

which provides an alternative expression to Eqs. (118), as: 



Cqq 


_ ' qq 


T 

1 -^qq 


C33 


V33 


rqq 




_ Vii- 


- Vlii 


C22 


V22- 


- VI22 



(121) 



V = ^Jl^Z^l , (122) 

t,33-i-3 
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that is equivalent to Eq. (114), and then admit real solutions provided in- 
equality (115) is satisfied. 

Finally, Eqs. (114) may be combined as: 

^11 _ Css'Pii ~ Cii'Pss . (123) 

^22 C337'22 — C22'P33 

where it can be seen that Eqs. (120) and (123) are changed one into the other, 
by replacing the terms, VasCqq/Cas^ with the terms, vlgg, and vice versa. The 
above results may be reduced to a single statement. 

Theorem 1. Given a R3 fluid, the relation: 

Xii Vii — Yii 



X22 V22 — Y22 
C33 



X,, = ^^33 , vigg ; g= 1,2 



C33 



yqq-vlqq , 7^7^33 ; 5=1,2 



is symmetric with respect to Xgg and Ygg, being the former tensor related 
to anisotropic random velocity distribution, and the latter one to systematic 
rotation around the axis, x^, or vice versa. 

In the special case of axisymmetric configurations, the dimensionless fac- 
tors appearing in the expression of the self potential-energy tensor, Eqs. (110), 
and the moment of inertia, Eqs. (Ill), do coincide with regard to equatorial 
axes, Vu = V22 andXn = X22, respectively, which necessarily imply di = (,22-1 
owing to Eq. (123). 

In the general case of triaxial configurations, the contrary holds, Vu 7^ 
V22 and Xii 7^ 222, then the equality, Cii — C22, via Eq. (120), implies the 
validity of the relation: 

„=^li:i^; (124) 

if otherwise, the random velocity distribution along the equatorial plane'' is 
anisotropic i.e. Cii 7^ C22- The related degeneracy can be removed using an 
additional condition, as it will be shown in the next section. 
The above results may be reduced to a single statement. 

^Throughout this paper, "along the equatorial plane" has to be intended as "along any 
direction parallel to the equatorial plane" . 
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Theorem 2. Isotropic random velocity distribution along the equatorial 
plane, Cii = C22, makes a necessary condition for R3 fluids to he symmetric 
with, respect to the rotation axis, X3. 

4 Imaginary rotation 

A unified theory of systematic and random motions is allowed, taking into 
consideration imaginary rotation. It has been shown above that Eq. (114), or 
equivalently one among (118), (122), admits real solutions provided inequal- 
ity (115) is satisfied. If otherwise, the rotation parameter, v, has necessarily 
to be negative, which implies, via Eq. (113), an imaginary figure rotation, iQ, 
where i is the imaginary unit. Accordingly, the related centrifugal potential 
takes the general expression: 



where Sgn is the sign function, Sgn(=p|a;|) = ^1, x ^ 0. The centrifugal 
force, dT/dw, is positive or negative according if real or imaginary rotation 
occurs, respectively. Then the net effect of real rotation is fiattening, while 
the net effect of imaginary rotation is elongation, with respect to the rotation 
axis (Caimmi 1996b; C06). 

To get further insight, let us particularize Eq. (114) to the special case of 
null rotation {v — 0). The result is: 



where the right-hand side, via Eqs. (31) and (110), depends on the mass 
distribution only. Accordingly, the net effect of positive (Cgg/Css > 0) or 
negative {(gg/Css < 0) random motion excess along the equatorial plane is 
fiattening (Vqq > V33) or elongation {Vqq < V33), respectively. In what 
follows, it shall be intended that random motion excess is related to the 
equatorial plane. 




2 2,2 
w = X1+X2 ; 

(125) 



Cqq_ _ 'Pqq 
C33 ^33 



V = 



(126) 
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4.1 Random motion excess and rotation 



In the limit of isotropic random velocity distribution, = ^22 = C33 — C/3, 
Eqs. (118) and (122) reduce to: 

Viso = ; q = l,2 ; (127) 

-Lqq 

v^so = ; (128) 

where the index, iso, means isotropic random velocity distribution. 
Accordingly, Eqs. (118) and (122) may be expressed as: 

V = ^^iso - ^^ani ; (129) 
^^ani=f^-ll^ ; q=l,2 ; (130) 



^iso ^ani j 










Iqq 




V33 







v,ni = h^ - 3 ; (131) 

\S,33 / M 

where Vani > for oblate- like configurations, Cgg/Css ^ Ij ^ani < for prolate- 
like configurations, C99/C33 ^ 1; index, ani, means contribution from 
random motion excess which, in general, makes an anisotropic random ve- 
locity distribution. Accordingly, positive or negative random motion excess 
is related to real or imaginary rotation, respectively. 
Let us rewrite Eq. (129) as: 

^^iso = V + ^^ani ; (132) 

which, owing to Eq. (113), is equivalent to: 

Ul^lf + Sgn(J--3^Ul, ; (133) 

where positive and negative Sgn vahies correspond to real and imaginary 
rotation, respectively. Then the effect of random motion excess on the shape 
of the system, is virtually indistinguishible from the effect of additional figure 
rotation. The above results may be reduced to a single statement. 
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Theorem 3. Given a R3 fluid, the effect of (positive or negative) random 
motion excess is equivalent to an additional (real or imaginary) figure rota- 
tion, Sgn((^/(^33 — 3)^2^^^^^, with regard to an adjoint configuration where the 
random velocity distribution is isotropic. 

Accordingly, a R3 fluid with assigned systematic rotation and random veloc- 
ity distribution, with regard to the shape, is virtually indistinguishable from 
an adjoint configuration of equal density profile, isotropic random velocity 
distribution, and figure rotation deduced from Eq. (133). 



4.2 Axisymmetric and triaxial configurations 

The combination of alternative expressions of the rotation parameter, Viso, 
defined by Eqs. (127), yields: 

X11P22-X22P11 = p33(Xii -X22) ; (134) 

which, for axisymmetric configurations, Tu — I22, Vu — 7^22, reduces to an 
indeterminate form, = 0. 

The combination of alternative expressions of the rotation parameter, 
fani, defined by Eqs. (129), yields: 

2^1 1 C22 ~ 2^22^11 = C33 (2^11 ~ 2^22) ; (135) 

which, for isotropic random velocity distributions, reduces to an indeter- 
minate form, = 0. In addition, axisymmetric configurations (In = I22) 
necessarily imply isotropic random velocity distributions along the equatorial 
plane, (ii = (22- 

The combination of Eqs. (107e) and (135) yields: 

^^^^C^^^^C^s^^^^^M . q^i^2 ; p = 2,l ; (136) 
which, for axisymmetric configurations (In —I22 =^3/2) reduces to : 

- ^(C - C33) = ; g = 1,2 ; (137) 

and the special case, C33 = C/3, reads Cii = C22 = C/3- 
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The limiting configuration, Qq = 0, via Eqs. (134), necessarily implies 
2pp < 2gq, owing to (pp > 0, g = 1, 2, p = 2, 1, and Eq. (136) reduces to: 

CIqq-C3s{2I,,-Ipp)^0 ; g = l,2 ; p = 2, 1 ; (138) 

which, owing to Eq. (107e), is equivalent to: 

^^(33^^2(33^ ; (139a) 

-^qq C33 C33 

Cqq = Q ; q=l,2 ; p = 2,l ; (139b) 

where Tqq/Tpp > 1 implies C33 > Cpp C33 ^ C/2- The above results may 
be reduced to the following statements. 

Theorem 4. Given a R3 Ruid, the anisotropy parameters along the equa- 
torial plane, Qq, q — 1,2, depend on the diagonal components of the dimen- 
sionless moment of inertia tensor, Iqq, q = 1,2, and the related expressions 
coincide, (u = C22, in the limit of axisymmetric conhgurations. In = l22- 

Theorem 5. Given a R3 fluid, a necessary and sufficient condition for isotropic 
random velocity distribution is that the anisotropy parameter along the ro- 
tation axis attains the value, C33 = 1/3. 

Theorem 6. Given a sequence of R3 fluids, the ending point occurs when 
the third diagonal component of the dimensionless self potential-energy ten- 
sor is null, V33 = 0, and/or the generalized anisotropy parameter related to 
the major equatorial axis is null, Cii = 0, which is equivalent to X22/X11 = 
(2 — C/C33)^^^- The related value of the rotation parameter is v = Vqq/Tqq, 
q — 1,2, independent of anisotropy parameters. The special case of dynami- 
cal (or hydrostatic) equilibrium, ( = 1, implies centrifugal support along the 
major equatorial axis, provided — 0. 

Accordingly, with regard to R3 fluids, the anisotropy parameters along the 
equatorial plane, Cii and C22, cannot be arbitrarily assigned, but depend on 

the dimensionless moment of inertia tensor diagonal components. In and 222, 
conform to Eqs. (136). On the other hand, the knowledge of the dimensionless 
moment of inertia tensor components, Xn and X22, the dimensionless self 
potential-energy tensor components, Vn, V22 and V33, together with the 
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rotation parameter, v, allows the determination of the rotation parameter, 
i^ani, via Eqs. (127), (128), (129), and then the ratios, Cgg/Css, Css/C, via 
Eqs. (130), (131), respectively, or the anisotropy parameter along the rotation 
axis, ^33, provided the virial index, (, defined by Eq. (107e), is assigned. 

In conclusion, with regard to R3 fluids deflned by assigned dimension- 
less moment of inertia tensor components, Xn and X22, dimensionless self 
potential-energy tensor components, Vu, V22 and V33, rotation parameter, 
V, and virial index, C, the anisotropy parameters, Cii, C22, CsS) cannot arbi- 
trarily be fixed, but must be determined as shown above. 

4.3 Sequences of virial equilibrium configurations 

With regard to R3 fiuids, it has been shown above that adjoints configura- 
tions are characterized by (i) centrifugal potential, 7iso(a;i, X2, X3) = T(xi, X2, 2^3)+ 

%ni{Xl,X2, X3), or ^l^{Xi, X2, X3) = n'^{Xi,X2, Xg) + Sgn(C/C33-3)^^am (Xi , 0:2, X3), 

Eq. (133); and (ii) isotropic random velocity distribution. Owing to Theorem 
3, a sequence of R3 fluids coincides with the sequence of adjoints conflgu- 
rations. Given a R3 fluid with fixed components of the dimensionless self 
potential-energy tensor, Vu, V22, V33, dimensionless moment of inertia ten- 
sor, Xii, X22, rotation parameter, v, and virial index, (, the anisotropy param- 
eters, ^11, ^22! C33) ai's determined via Eqs. (127)-(131) and (136). Negative 
values of the rotation parameter, Viso, extend the sequence of axisymmetric 
configurations to imaginary rotation i.e. prolate configurations where the 
major axis coincides with the rotation axis. 

Once more owing to Theorem 3, the bifurcation point of a sequence of R3 
fluids coincides with its counterpart along the sequence of adjoint conflgura- 
tions. Aiming to flnd a necessary condition for the occurrence of a bifurcation 
point, let us equalize the alternative expressions of Eqs. (118). The result is: 



where the anisotropy parameter ratios, CqqK?>?>i 1 = 1)2, may be deduced 
from Eqs. (136). After some algebra, Eq. (140) reads: 



Vn 
In 




(140) 



2^117^22 ~ 2^22^11 



^33 



(141) 
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and the occurrence of a bifurcation point has necessarily to satisfy the rela- 
tion: 

TnP22 - T22Vn ^ . ^^^^^ 

£21-^1 111 - -L22 

where €21 = 02/^1 is the ratio of two generic equatorial (perpendicular) radii, 
and £21 1 implies X22 — > Xn, P22 ^ ^ii- Then Eq. (142) is a necessary 
condition for the existence of a bifurcation point, as it selects the sole axisym- 
metric configuration which satisfies Eq. (141), regardless from the values of 
the anisotropy parameters, Cr-r, '"=1,2, 3. The above results may be reduced 
to a single statement. 

Theorem 7. Given a sequence of R3 fluids, a necessary condition for the 
existence of a bifurcation point from axisymmetric to triaxial configurations, 
is independent of the anisotropy parameters, („, ^ = 1,2,3, and coincides 
with its counterpart related to the sequence of adjoint configurations with 
isotropic random velocity distribution. 

Sequences of R3 fluids with assigned density profiles, can be deduced from 
the knowledge of the rotation parameters, i^iso(e3i) and t'ani(e3i, Csa), as func- 
tions of the meridional axis ratio, €31 = a^/ai, and the effective anisotropy 
parameter, (^33, as represented in Figs. 2 and 3, respectively^. 

A hypothetical sequence of axisymmetric R3 fiuids, extended to the case 
of imaginary rotation, is shown in Fig. 2. Hypothetical dependences of the 
rotation parameter, t;ani, on the meridional axis ratio, £31, and the effective 
anisotropy parameter, ^33 = C33/C (labelled on each curve), are shown in 
Fig. 3. With regard to a fixed effective anisotropy parameter, the related 
sequence starts from a nonrotating configuration, v = i.e. v^m = Viso, and 
ends at a configuration where 631 = and/or (n = 0. A hypothetical locus 
of the ending points related to the latter alternative, is represented as a long- 
dashed curve. The locus of nonrotating configurations (short-dashed lines) 
coincides with the curves represented in Fig. 2. No sequence can be continued 
on the right, as imaginary rotation i.e. larger ^33 would be needed and a 
different sequence should be considered. The effect of positive (^33 < 1/3) or 
negative (C33 > 1/3) random motion excess is equivalent to an additional real 

^Strictly speaking, the curves plotted in Figs. 2, 3, and 4 are related to the special case 
of homeoidally striated Jacobi ellipsoids (C06), but can be taken as illustrative for R3 
fluids. 
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Figure 2: A hypothetical sequence of axisymmetric R3 fluids, from the start- 
ing point (square) to the bifurcation point (triangle), and related triaxial 
R3 fluids, from the starting point (triangle) to the bifurcation point (Greek 
cross) , extended to the case of imaginary rotation (negative values of the rota- 
tion parameter, v). In any case, the random velocity distribution is isotropic. 
Both sequences are continued in the region of instability. The horizontal line, 
Visa = 0, is the locus of non rotating and/or zero volume configurations. The 
vertical line, £31 = 1, is the locus of round(ai = 02 = 03) configurations. The 
above mentioned lines divide the non negative semiplane, £31 > 0, into four 
regions (from top left in clockwise sense): A - oblate-like shapes with real ro- 
tation; B - prolate-like shapes with real rotation; C - prolate-like shapes with 
imaginary rotation; D - oblatc-Iikc shapes with imaginary rotation. Regions 
B and D are forbidden to sequences 5rR3 fluids. 



_ 1 1 1 1 1 1 1 1 




- 1 1 1 1 1 1 1 1 


\ ; 



0.0 0.5 1.0 1.5 

^31 



Figure 3: A hypothethical dependence of the rotation parameter, -Uani, related 
to random motion excess, on the meridional axis ratio, eai, with regard to 
R3 fluids. Each curve is labelled by the value of the effective anisotropy 
parameter, ^33 = C33/C- The horizontal non negative semiaxis, £31 > 0, Vani = 
0, is the locus of configurations with isotropic random velocity distribution, 
C33 = 1/3- The vertical straight line, 631 = 1, is the locus of round (ai = 
0'2 — 0,3) configurations. The generic sequence starts from a non rotating 
configuration (short-dashed lines) and ends at a configuration where either 
£31 = and/or (u — (long-dashed line). The regions above the upper short- 
dashed curve and below the long-dashed curve, respectively, are forbidden to 
R3 fluids. The non negative vertical semiaxis, v^ni > 0, 631 = 0, corresponds 
to fiat (C33 = 0) configurations with no figure rotation. The curves are 
symmetric with respect to the horizoS^l axis, until the limiting curve, (33 = 
C = 1, is reached. The limiting configuration where Cn = 0, is marked by 
a square: no configuration in virial equilibrium is allowed on the left, as it 
would imply < 0. No configuration is allowed on the right of the starting 
point, as it would imply v < 0. 



or imaginary rotation, respectively. The horizontal non negative semiaxis, 
esi > 0, Dani = 0, is the locus of configurations with isotropic random velocity 
distribution, C33 = 1/3. The vertical straight line, £31 = 1, is the locus of 
round ( ) configurations. 

With regard to real rotation, the ending configuration (631 = and/or 
(11 = 0) is marked by a square. Configurations on the left arc forbidden, as 
they would imply negative random energy tensor component, (-Erdm)!! < 0, to 
satisfy the virial equation (107a), which demands imaginary rotation around 
major equatorial axis i.e. systematic motions other than rotation around the 
minor axis. 

With regard to imaginary rotation, no ending point occurs and the system 
is allowed to attain negative infinite rotation parameter, Viso —00, and 
infinite meridional axis ratio, £31 — > +00. The related configuration is either 
a spindle (oi = a2 = 0) or an infinitely high cylinder (oi = 02 < 03 — > +00). 

Further inspection of Fig. 3 shows additional features, namely: (i) null left 
first derivatives on each sequence at bifurcation point (not marked therein for 
sake of clarity), and (ii) occurrence of symmetric sequences with respect to 
the horizontal axis (including also forbidden configurations). For additional 
considerations on (i), see C06 (nothing changes with respect to the special 
case investigated therein). The above results may be reduced to a single 
statement. 

Theorem 8. Given a sequence of R3 fluids, the contribution of random mo- 
tion excess, Vam, to the rotation parameter, Viso, has a null left first derivative 
at the bifurcation point, [d-Uani/ de3i](e3j)-. = 0. 

The occurrence of symmetric sequences with respect to the horizontal 
axis, is deduced from Eq. (131) using the condition: 

((3^3)-' - 3 = -(C3-3)-' + 3 ; (143) 

where (^;^ = C33/C is related to curves characterized by negative (^33 = C33 ^ 
1/3) and positive (C33 = C33 ^ 1/3) values, respectively, of the rotation pa- 
rameter, t^ani, see Fig. 3. Couples of symmetric sequences (including forbidden 
configurations) start from (C33,C33) = (1/3, 1/3), where each curve coincides 
with the other, and end at (1,1/5). Sequences in the range, < C33 < 1/5, 
have no symmetric counterpart. 

Let a point, P[e3i, 'Uiso]) be fixed on a sequence of axisymmetric R3 fiuids, 
and the dimensionless moment of inertia tensor components. In and X22, 
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be determined by use of Eqs. (127) and (128). Let a point, P'(e3i, t;ani)j be 
fixed on the plane, (0 €31 t'ani), and the effective anisotropy parameters, 
Cii) C22, CsS) be determined by use of Eqs. (130), (131), (136). Finally, let the 
rotation parameter, v, be determined by use of Eq. (129). Following the above 
procedure, sequences of R3 fluids may be generated. For assigned density 
profiles and systematic rotation velocity fields, there are three independent 
parameters, which may be chosen as two axis ratios, £21, £31, and one effective 
anisotropy parameter, (^33. 

In the (Oesiv) plane (Fig. 3), each sequence starts from the intersection 
between the curves, v = [';j(e3i)]iso (dashed), v = ["^(£31; C33)]ani (full), and 
ends at the intersection between the curves, v = v{esi; (u = 0) (long-dashed), 

V = b(e3i;C33)]am (Ml). 

Hypothetical dependences of the rotation parameter, v, on the meridional 
axis ratio, £31, and the effective anisotropy parameter, ^33 = C33/C (same cases 
as in Fig. 3), are shown in Fig. 4. A hypothetical locus of the ending points 
related to (u = is represented as a long-dashed curve, corresponding to 
1/2 < (^33 < 1. The continuation of a generic sequence on the left of the long- 
dashed curve would imply < or £31 < 0. The ending point of sequences, 
related to < C33 < 1/2, coincides with the origin. The initial configuration, 
related to ^33 = 1, corresponds to = Oi = 02 < 03 or Oi = 02 < 03 — > +00, 
which is equivalent to £31 +00. 

4.4 A special case: homeoidally striated Jacobi ellip- 



Homeoidally striated Jacobi ellipsoids are a special case of R3 fluids, for which 
the results are already known (CM05; C06). Then the particularization of 
the current theory to homeoidally striated Jacobi ellipsoids makes an useful 
check. In the case under discussion, Eqs. (110) and (112) reduce to (e.g.. 



soids 



CMOS; C06): 




(144) 



rr 



r = 



1,2,3 ; 




r = l,2,3 ; 



(146) 



(145) 
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Figure 4: A hypothethical dependence of the rotation parameter, -u, related to 
systematic rotation, on the meridional axis ratio, eai, with regard to R3 fluids. 
Each curve is labelled by the value of the effective anisotropy parameter, 
C33 — C33/C) except the lower two, where C33 = 2/9,1/5, respectively. The 
dotted vertical straight lines denote a hypothetical bifurcation point (left) 
and the round ( ) configuration. The generic sequence starts from 

a non rotating configuration on the horisontal axis and ends at a configuration 
where £31 = or = (long-dashed fine), denoted by a square. The 
continuation on the left of the ending point, where Cii < 0, 1/2 < (^33 < 1, 
is also shown. The initial configuration, related to C33 = 1, corresponds to 
= ai = a2 < 03 or ai = a2 < 03 — > +00, which is equivalent to 631 — > +00. 
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where Ugci is a profile factor i.e. depending only on the radial density pro- 
file, and Ar, r = 1,2, 3, are shape factors i.e. depending on the axis ratios 
only, and a^, r = 1,2,3, are the semiaxes of the ellipsoidal boundary. The 
dimensionless density profile may be represented as: 

p = po/(0 ; /(i) = i ; Po = p(i) ; (i47a) 
e = o<e<5; H = -; (i47b) 

ro ro 

where po, Tq, are a scaling density and a scaling radius, respectively, with 
respect to a reference isopycnic surface, while S = E{R,9,(j)), and R are 
related to the truncation isopycnic surface. 

The mass, M, and the moment of inertia tensor, Ipq, are (e.g., CMOS): 

M = z/„,asMo ; (148) 
Ipq = 5pqVi^^Mal ; (149) 

where Mq is the mass of a homogeneous ellipsoid with same density and 
boundary as the reference isopycnic surface, and z^masj ^'im-j ai'c profile factors. 
The combination of Eqs. (102), (111), (112), (147b), and (149) yields: 

= 2:11 + X22 = Vinr -621631 (1 + 4) ; (151) 



which allow the particularization of the general results related to R3 fluids, 
to homeoidally striated Jacobi ellipsoids. 

Using Eqs. (118), (145), (150), and (151), and performing some algebra, 
the rotation parameter, v, takes the expression: 

^^2_u^ C,,A,-Cqqe%A, q^i^2 ■ (152) 

3 Vinr C33 

let us deflne a normalized rotation parameter, as: 

= V ; (153) 
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accordingly, Eq. (152) reduces to: 



C33 



Vn = 7 — ; g = l,2 ; (154) 



which, in spite of a different definition of the rotation parameter, v, coincides 
with a previously known result (COG) i.e. — (fN)co6- 

Using Eqs. (119), (146), (150), and performing some algebra, the anisotropy 
parameter ratio, Cgg/CsS) takes the expression: 

which, keeping in mind a different definition of the rotation parameter, after 
some algebra can be shown to coincide with a previously known result (COG). 
In addition, Eqs. (155) disclose that: 

Cll ^1 - ^^N 

which, owing to Eqs. (155), has necessarily to coincide with a previously 
known result (COG). 

Using Eqs. (122), (14G), (151), and performing some algebra, the alterna- 
tive expression of the rotation parameter, v, takes the form: 

U + e2ljC33 

which, in spite of a different definition of the rotation parameter, v, coincides 
with a previously known result (COG) i.e. = (i'n)co6- 

Using Eqs. (123), (145), (151), and performing some algebra, the di- 
mensinless moment of inertia tensor component ratio, 222/2^11, takes the 
expression: 

2^22 _ 2 _ C33e21^2 — C22e31^3 roN 

^ - ^21 - — ^^24 ; (^5Sj 

Ml 1,33^1 — ^lie31^3 

which, after additional algebra, can be shown to coincide with a previously 
known result (COG; the counterpart of X22/2^ii is 7^22/^11 therein). 
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Using Eqs. (124), (145), (150), (153), and performing some algebra, the 
rotation parameter, v, related to isotropic random velocity distribution along 
the equatorial plane, Cii = C22, takes the expression: 

= — — 2 — ; (1^9) 

^ ~ ^21 

which, after additional algebra, can be shown to coincide with a previously 
known result (COG). 

The above results hold, in particular, for isotropic random velocity dis- 
tributions {C,ii = C22 = C33 = C/3), which implies that the expressions of the 
rotation parameters, Viso and i^ani, via Eq. (132), must necessarily coincide 
with their previously known counterparts (COG). 

Using Eqs. (13G), (150), (151), and performing some algebra, the anisotropy 
parameters along the equatorial plane, Qg, q — 1,2, take the expression: 



Ce^l - C33(26^i - 6^ 
1 + ^21 



Qq - ^ 7 2 — - ; ? = i,2 ; (iGo) 



which coincide with previously known results (COG), in particular for axisym- 
metric configurations (e2i = 1)- 

Using Eqs. (139) and (158), the condition for the occurrence of the ending 
configuration, related to Cii = 0, takes the expression: 



2 C33 — C22 2(^33 — ( 



6^1 = ^^^^ = ^"^ ; Cii = ; (IGl) 

<533 S,33 

which coincides with a previously known result (COG). 

Using Eqs. (141), (145), (150), and performing some algebra, a general 
relation between dimensionlcss self potential-energy tensor components and 
dimensionless moment of inertia tensor components, takes the expression: 



6^1(^2-^1) 



= ; (162) 



(1 - 4: 

which coincides with a previously known result (Caimmi 199Ga; COG). Ac- 
cordingly, a necessary condition for the occurrence of a bifurcation point, 
Eqs. (142), reduces to: 

li- ^^^^S^ = 4^3 ; (163) 

£21— >1 i — 

which coincides with a previously known result (Caimmi 199Ga; COG). 
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5 Conclusion 



The current paper has been aimed in getting more insight on three main 

points concerning large-scale astrophysical systems, namely: (i) formulation 
of tensor virial equations from the standpoint of analytical mechanics; (ii) 
investigation on the role of systematic and random motions with respect 
to equilibrium configurations; (iii) extent to which systematic and random 
motions are equivalent in flattening or elongating the shape of a mass distri- 
bution. 

The tensor virial equations have been formulated regardless from the 
nature of the system and its constituents, by generalizing and extending a 
procedure used for the scalar virial equations, in presence of discrete subunits 
(Landau & Lifchitz 1966, Chap. II, §10). In particular, the self potential- 
energy tensor has been shown to be symmetric with respect to the indices, 
{Epot)pq = {Epot)qp- Then the results have been extended to continuous mass 
distributions. 

The role of systematic and random motions in coUisionless, ideal, self- 
gravitating fluids, has been analysed in detail including radial and tangential 
velocity dispersion on the equatorial plane, and the related mean angular 
velocity, fi, has been conceived figure rotation. 

R3 fluids have been defined as ideal, self-gravitating fiuids in virial equi- 
librium, with systematic rotation around a principal axis of inertia. The 
related virial equations have been written in terms of the moment of iner- 
tia tensor, Ipg, the self potential-energy tensor, {Epot)pq, and the generalized 
anisotropy tensor, (pg (CMOS; C06). Additional effort has been devoted to 
the investigation of the properties of axisymmetric and triaxial configura- 
tions. 

A unified theory of systematic and random motions has been developed 
for R3 fluids, taking into consideration imaginary rotation (Caimmi 1996b; 
C06). The effect of random motion excess has been shown to be equivalent 
to an additional real or imaginary rotation, inducing flattening (along the 
equatorial plane) or elongation (along the rotation axis), respectively. Then 
it has been realized that a R3 fluid always admits an adjoint conflguration 
with isotropic random velocity distribution. 

In addition, further constraints have been established on the amount of 
random velocity anisotropy along the principal axes, for triaxial conflgu- 
rations. A necessary condition has been formulated for the occurrence of 
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bifurcation points from axisymmetric to triaxial configurations in virial equi- 
librium, which is independent of the anisotropy parameters. 

The particularization to the special case of homeoidally striated Jacobi 
ellipsoid has been made, and some previously known results (C06) have been 
reproduced. 
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Appendix 

A. Tensor potentials and potential-energy tensors 

For reasons of simplicity (integrals are easier than summations to be cal- 
culated), let us take into consideration a continuous distribution of matter, 
where volume elements, A^*, interact each with the other according to their 



43 



charge density, = A^^/AS*, being A^^^ the charge within AS*. It is in- 
tended that the results found in this section can be extended to discrete mass 
distributions, using summations instead of integrals. Let (0 Xi X2 X3) be a 
reference frame where the origin coincides with the centre of mass, and the 
coordinate axes with the principal axes of inertia. 

The effect of the interaction on an infinitesimal volume element, d^S* = 
dxi da;2 dxs, due to a charge distribution of density, p^{xi,X2,X3), is deter- 
mined by the potential: 



.5=1 



(164) 



where the constants, x ^^^^ G^,, specify the nature and the intensity of the 
interaction, respectively, and S is the volume filled by the system. 

The first derivatives of the potential with respect to the coordinates, are: 



dV 

dx. 



P-^{Xy^X2-,X2){Xs 

H l-x/2 



■ 3 
.s=l 



' \2 



s = 1,2,3 



(165) 



it can be seen that the functions of the coordinates, V and XpdV /dxq, 1 < 
P<3,l<g<3, are homogeneous functions of degree x, and the Euler 
theorem holds (e.g., LL66, Chap. IV, §10). 

With regard to a selected infinitesimal volume element, the potential 
may be thought of as the tidal energy due to the whole charge distribution, 
related to the point under consideration, with the unit charge placed therein. 
Associated with the potential, defined by Eqs. (164) and (165), is the self 
potential energy: 



Epot = ~ J^Pxi^l^^2,X3)V{Xi,X2,X3)d^S ] 



(166) 



the self potential energy may be thought of as the tidal energy due to the 
whole charge distribution, related to all the infinitesimal volume elements, 
provided any pair is counted only once. 

The coincidence of Eqs. (166) and (32) may be verified along the following 
steps: (i) write the alternative expressions of the potential energy in explicit 
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form, using Eqs. (164) and (165); (ii) express the explicit form of Eq. (32) 
as a sum of two halves; (iii) with regard to one half, replace the variables 
of integration, {xi,X2,X3) <-> {x[,X2,x'^), keeping in mind that the integrals 
are left unchanged; (iv) sum the resulting two halves and compare with the 
explicit form of Eq. (166). 

To get more information on the charge distribution, let us define the 
tensor potential: 



Vpq{Xi,X2,X3) =G^j p^{x\,x'2,x'3) 



[Xp Xp) [Xq 



x„ 



.s=l 



l-x/2 



(167) 



and the self potential-energy tensor: 

(-Epot)pg ^ Pxi^l^ ^2, X3)Vpq{Xi, X2, X3) S' 

the above mentioned tensors are manifestly symmetric: 

Vpq{Xi,X2, Xs) = Vqp{Xi,X2, X3) ; 
{Epot)pq — {Epot) qp 't 

and the related traces equal their scalar counterparts: 

3 

^ Vss{xi,X2, X3) = V{xi,X2, X3) ; 
s=l 

3 

y^XEpot) ss = Epot ; 



(168) 



(169) 
(170) 



(171) 
(172) 



conform to Eqs. (164), (167), and (166), (168), respectively. 

The coincidence of Eqs. (168) and (31) may be verified along the follow- 
ing steps: write the alternative expressions of the potential-energy tensor in 
explicit form, using Eqs. (167) and (165); (ii) express the explicit form of 
Eq. (31) as a sum of two halves; (iii) with regard to one half, replace the 
variables of integration, {xi,X2,X3) {x[,X2,x'^), keeping in mind that the 
integrals are left unchanged; (iv) sum the resulting two halves and compare 
with the explicit form of Eq. (168). 
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In the special case of gravitational interaction, % = — 1, = G (constant 
of gravitation), Px = P (mass density), the potential and the potential energy, 
both in scalar and in tensor form, attain their usual expressions known in 
literature [C69, Chap. 2, § 10; see also therein a proof for the equivalence 
of Eqs. (31) and (168) - also illustrative for the equivalence of Eqs. (32) and 
(166) - where the steps outlined above are followed in a reversed order]. 



B. An alternative expression of the generalized anisotropy 
parameters 

Let us rewrite Eqs. (107a) and (107b) in a more compact notation, as: 

(1 - 5^r)Irr^^ + Ma^, + (^pot)r-r = Ma^, - MCrr(J^ ; T = 1, 2, 3 ; (173) 
and the combination of Eqs. (105a) and (173) yields: 

U„^M{al-CrrCT^) ; r = l,2,3 ; (174) 

from which the following expression of the generalized anisotropy parameters 
is derived: 

and the related trace, owing to Eqs. (15), (107c), and (107e) reads: 

where C, exceeds unity for a moment of inertia with respect to the centre of 
mass decreasing in time, and vice versa. 
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